Multigrid Method for H(div) in Three Dimensions∗

We are concerned with the design and analysis of a multigrid algorithm for H(div; Ω)–elliptic linear variational problems. The discretization is based on H(div; Ω)–conforming Raviart–Thomas elements. A thorough examination of the relevant bilinear form reveals that a separate treatment of vector fields in the kernel of the divergence operator and its complement is paramount. We exploit the representation of discrete solenoidal vector fields as curls of finite element functions in so-called Nédélec spaces. It turns out that a combined nodal multilevel decomposition of both the Raviart–Thomas and Nédélec finite element spaces provides the foundation for a viable multigrid method. Its Gauß–Seidel smoother involves an extra stage where solenoidal error components are tackled. By means of elaborate duality techniques we can show the asymptotic optimality in the case of uniform refinement. Numerical experiments confirm that the typical multigrid efficiency is actually achieved for model problems.